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ABSTRACT 






C^^ ' Context. Interstellar magnetic fields influence all stages of the process of star formation, from the collapse of molecular cloud cores 

to the formation and evolution of circumstellar disks and protostellar jets. This requires us to have a full understanding of the physical 
properties of magnetized plasmas of different degrees of ionization for a wide range of densities and temperatures. 
Aims. We derive general equations governing the magneto-hydrodynamic evolution of a three-fluid medium of arbitrary ionization, 
also including the possibility of charged dust grains as the main charge carriers. In a companion paper (Pinto & Galli 2007), we 
complement this analysis computing accurate expressions of the collisional coupling coefficients for a variety of gas mixtures relevant 
Q , for the process of star formation. 

5^ . Methods. Over spatial and temporal scales larger than the so-called large-scale plasma limit and the collision-dominated plasma limit, 

"^ ' and for non-relativistic ffuid speeds, the electric field, the electric current and the ion-neutral drift have their instantaneous values de- 

rrt ' termined by the evolution of the magnetic field, which obeys an advection-diffusion equation. The validity of the approximations 

made is discussed critically. 

Results. We derive the general expressions for the resistivities, the diffusion time scales and the heating rates in a three-fluid medium 
and we use them to estimate the evolution of the magnetic field in molecular clouds and protostellar jets. Collisions between charged 
^ ' particles significantly increase the value of the Ohmic resistivity during the process of cloud collapse, affecting in particular the 

^O , decoupling of matter and magnetic field and enhancing the rate of energy dissipation. The Hall resistivity can take larger values 

04 ■ than previously found when the negative charge is mostly carried by dust grains. In weakly- or mildy-ionized protostellar jets, am- 

00 ' bipolar diffusion is found to occur on a time scale comparable to the dynamical time scale, limiting the validity of steady-state and 

04 ' nondissipative models to study the jet's structure. 

"^ ' Key words. MHD;Plasmas;ISM:magnetic fields;ISM:clouds;ISM:jets and outflows 
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^ . and several charged atomic and molecular species. About 2% of ^q^ cm-^ if a large number of very small grains is present (Nishi, 

■ - - the mass of the ISM ism the form of solid particles (dust grains) Nakano & Umebayashi 1991), or during the collapse of molec- 

with sizes ranging from ~ 10 A to- 1 jum that may carry electric ular clouds when the neutral density becomes larger than ~ 

charge (Spitzer 1978, Nakano & Umebayashi 1980, 1986). The lo^ cm"^ (Nakano et al. 2002). Second, many of the commonly 

transformation of a piece of ISM into a protostar involves vari- adopted expressions for the conductivity of the ISM in collaps- 

ations of several orders of magnitude in density, temperature, ing clouds and protostellar jets are derived neglecting collisions 

and ionization fraction, resulting in a variety of interconnected between charged particles (e.g., Wardle & Ng 1999; Desch & 

environments, such as cold dense clouds, cicumstellar disks, and Mouschovias 2001 ; Nakano, Nishi & Umebayashi 2002; Tassis 

protostellar jets. Magnetic fields affect the motion of the charged & Mouschovias 2005). This so-called "multifluid approach" can 

component of the gas, which in turn transfers the effects of the be used for an arbitrary number of charged species of arbitrary 

Lorentz force to the neutrals by elastic collisions. mass, but is valid only for low-ionization degrees and cannot be 

While several authors have focused on specific ionization generalized to include collisions between charged particles in a 

regimes or chemical compositions, we lack a rigorous formula- self-consistent way (see e.g., Kelley 1989). Benilov (1996, 1997) 

lion of the magneto-hydrodynamic (MHD) equations in general derived the multifluid equations for a nonmagnetized plasma in 

valid for the variety of situations occurring in the star formation detail, including accurate expressions for the rate of momentum 

process. First, the standard expressions of the plasma conduc- transfer by elastic and inelastic collisions, 

tivities are generally derived neglecting the electron's mass with Following a different approach, in this paper we derive the 

respect to the ion and neutral's masses (e.g., Bragmski 1965, gg^gj^i equation for the evolution of the magnetic field in a gas 

of arbitrary ionization degree made of positively- and negatively- 

Send offprint requests to: C. Pinto charged particles and neutrals. We also obtain the general ex- 



1. Introduction Mitchner & Kruger 1973), and are not applicable to condi- 

tions where the carriers of negative charge are dust grains 
The interstellar medium (ISM) is a partially-ionized gas contain- rather than electrons. For example, charged dust grains domi- 
ing substantial fractions of neutral particles, mosdy H and H2, nate over free electrons in clouds cores of neutral density ~ lO'* 



Pinto et al.: Three-fluid MHD equations 



pressions of the resistivity coefficients, the associated timescales, 
and the heat generation rates for a three-fluid plasma. In a com- 
panion paper (Pinto & GalH 2007, hereafter Paper II), we com- 
pute accurate numerical and analytical expressions for the colli- 
sional coupling coefficients that are required to calculate the re- 
sistivity coefficients of the plasma. As an illustration, we apply 
the formalism derived in this paper to study the evolution of the 
magnetic field in a molecular cloud and in a protostellar jet with 
the help of simple analytical models. The same formalism can 
be adopted to study other situations such as e.g., multicompo- 
nent radiatively-driven stellar winds (Krticka & Kubat 2001) and 
shock waves in plasmas of low-ionization degree (Mullan 1971, 
Draine 1980, Draine & McKee 1993; Guillet, Pineau des Forets 
& Jones 2007). 

In Sect. |2]we develop a rigorous three-fluid theory for plas- 
mas of arbitrary degree of ionization and particle masses, crit- 
ically discussing the approximations made in deriving the gov- 
erning equations; in Sect. |3]we obtain the complete advection- 
difiiision equation for the magnetic flux for an axially-symmetric 
system; in Sect. |4]we obtain the equation for the rate of energy 
generation due to friction forces between the fluid components; 
in Sect. |5] and in Sect. |6] we apply the three-fluid equations to 
typical conditions of molecular clouds and protostellar jets, re- 
spectively, evaluating the relevant evolutionary timescales and 
the heating rates; finally, in Sect. |7] we summarize our conclu- 
sions. 



2. Three-fluid description of a partially-ionized gas 

We consider a system composed by three fluids, namely 
positively-charged particles (subscript +), negatively-charged 
particles (-), and neutrals (n). The generic species s (with s - 
+, -, n) is characterized by particle mass m,, density ps, mean 
flow velocity u, and charge ^j, with q^ - -Z e and q+ - Z+e, 
where e is the electron charge. Charged particles may include 
electrons (e), ions (/) and positively- or negatively-charged dust 
grains ig^,g^)- We define the quantity 



Z+m^ 



(1) 



the mass ratio of negatively- and positively-charged species. 
No assumption is made on the magnitude of e: for a plasma 
of ions and electrons, e « 1; for a plasma of ions and 
negatively-charged grains, e » 1 ; for a plasma of negatively- 
and positively-charged grains, e a; 1 . 

In the following, and in Paper II, we specialize to the case 
where the fluid components are characterized by a maxwellian 
velocity distribution function with temperatures T^, that is not 
necessarily the same for all species. Necessary conditions for 
this assumption are that (/) the characteristic evolution timescale 
of the system must be much longer than the inverse of the colli- 
sion frequency between particles of a given species and (//) the 
drift velocity between particles of different species must be suffi- 
ciently small to produce a frictional heating rate smaller than the 
thermal energy density of each component divided by its self- 
relaxation time (Draine 1986). The assumption of a maxwellian 
velocity distribution is equivalent to assume that the stress tensor 
is diagonal and, therefore, to neglect viscosity and heat conduc- 
tion in the fluid (see Sect. 12. 11 ). 

The plasma is also electrically neutral. 



2-^^ = 0, 



(2) 



or, by Eq. ([T]), 
P- = ep+- 



(3) 



The assumption of charge neutrality is valid if the characteristic 
length scale of the system is much larger than the Debye length 
of the charged particles (see Sect. l2.8l and Paper II). 

2.7. Definitions 

The mean velocity of each species averaged over its maxwellian 
velocity distribution function fsi\s) is given by 



U, = {\}s 



/v,/j(Vj)(iv, 

JfsiVs) d\,, 



(4) 



where Vj is the microscopic velocity of the particles, and fs(\) is 
normalized as 



/ 



/,(v)flfv=l. 



We also define the mass fractions of the three fluids, 

hs — > 

p 



(5) 



(6) 



with Yjs^s - 1- The average density and velocity of the three- 
fluid system are 



P^YjP-. 



and 

U=-J]p,u„ 



(7) 



(8) 



respectively. 

It is convenient to identify the random component of the mi- 
croscopic velocity of each species with respect to the average 
fluid velocity U. Therefore, we define 



V.S = U + W.S, 

so that 

U.S = (U + W>.5 = U + (w).5. 



(9) 



(10) 



where { )s indicates the average performed on the distribution 
function of the particles of species s. 

Consistent with the assumption of a maxwellian velocity dis- 
tribution, the pressure tensor for each species Pj = Ps{wjWk)s is 
diagonal. 



Pf^Ps6jk, with P,^-Ps{w^}s- 

It is possible then to define the kinetic temperature T^ 

3 1 2 

-kBTs = -m,{w )s, 



(11) 



(12) 



where ks is the Boltzmann's constant, and the thermal speed a, 
for each species 



2 ^bTs 

a, = , 



(13) 



in terms of the mean random energy of each component. Notice 
that this definition of temperature is based on the microscopic 
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velocities of each species with respect to the average fluid, not 
with respect to the mean velocity of each species. 

With these definitions, the pressure of each component is 
then given by 



Ps = a^ps 



(14) 



Since the random components Wj of the velocity of the different 
species are defined with respect to the same frame of reference 
(the one in which the center of mass is at rest), the pressure con- 
tributions of each species can be consistently added together to 
form the total pressure of the average fluid. 



p^Y.'. 



^p, where a^ = y ^sC 



(15) 



We also define the total density and pressure of the charged 

species as 



Pc=P++P- 



P, = P^ + P^. 



(16) 



We stress that this formulation, commonly adopted in 
plasma physics, (see, e.g., Krall & Trivelpiece 1973, Boyd & 
Sanderson 1969), is different from the approach adopted in the 
astrophysical literature (e.g., Mouschovias 1991)where the ran- 
dom velocity of each species is defined as a deviation with re- 
spect to the mean velocity of the same species, v^ - u^-i-Wj. With 
the latter definition, the derivation of the fluid equations for each 
species is simplified by the fact that the average of the random 
components (w')., are zero. However, the pressure and the heat 
tensor defined in this way have a different meaning with respect 
to the corresponding quantities in our derivation. Although the 
partial pressures defined in the two approaches are nearly equiv- 
alent when Us ~ U, the definitions adopted in the two cases, 
however, differ from a conceptual point of view, because in our 
approach the random components of the velocity of each species 
are defined in the same reference frame, and the partial pressures 
can be added to form physically meaningful global quantities. 

By eq. ([3]), the electric current is 



—psUs = P-(U+ - U_) 

m. „ m. 



Following Schliiter (1950, 1951), we also define an "ambipolar 
current" J^ associated to the ion-neutral drift speed u+ - u„. 



1- 

Jd = Pe(u+ 



u„). 



(18) 



2.2. Equations of continuity and momentum 



Taking the zeroth and first order momenta of the Boltzmann 
equation for each species, the equations of continuity and mo- 
mentum read: 

dps 
dt 



^+V-(p,u,) = ^5, 



(19) 



dt 



(p,u,) + V ■ (p,<vv),) = -p.W 



(20) 



where 'V is the gravitational potential, E is the electric field, B 
is the magnetic field, F,.,/ is the friction force (per unit volume) 
exerted on particles of the fluid s by particles of the fluid s' by 
elastic collisions, S ss' and r 5^/ are the rate of change of density 
and momentum of the species s due to inelastic collisions (chem- 
ical reactions) with particles of species s'. 



2.3. Elastic coliisions 

The calculation of the friction force Fjs' accounting for elastic 
collisions between particles with a Maxwellian velocity distri- 
bution is reviewed in Paper II. The resulting expression is 

F.V.S' =a'jj'(Uj' -Uj), (21) 

where a,,,' is the so-called/ncf!o« coefficient, related to the mo- 
mentum transfer rate coefficient {crv}ss' by 

PsP.s 



m^ + m. 



-{o-v)s 



(22) 



For elastic collisions, Uss' - cts's, and therefore Fjy — -Fj-j. In 
general, {a-v)ss' is a function of the temperatures Tj and T,' and 
of the relative mean velocity |u.s' - Us| of the interacting species. 
Thus, Fjs' is a nonlinear function of the relative mean velocity. 
Accurate values of {crv)ss' for various colliding particles, and 
useful fitting formulae for a wide range of temperatures and drift 
velocities, are derived in Paper II. 

Another frequently used quantity is the collisional drag co- 
efficient 

_ {(rv)ss' ,,„. 

Jss' = , (23) 

m.5 + ms' 

which is independent on the particle's density. 



2.4. Inelastic collisions 

The expressions for the rate of change of momentum due to in- 
elastic collisions depend on the specific type of chemical reac- 
tion considered. For binary reactive collisions that produce one 
or two particles, Benilov (1996) has found that 



, inet 



a„, -(uy -Uj)-h5,,'U,, 



(24) 



where c'™' is a complex function of the masses and tempera- 
tures of the two reacting species, their relative mean velocity, and 
the reaction cross section (or the reaction rate). Approximated 
(1 ') expressions for a™' have been obtained by Draine (1986) for 
several specific reactions, and by Ciolek & Mouschovias (1993) 
for reactions involving dust grains. 

Within a three-fluid scheme like the one considered in this 
paper, the only possible reactions between two species of par- 
ticles are recombinations of singly-charged ions with electrons, 
resulting in the production of neutral particles (i + e — » n) and 
ionization reactions (« — > / + e) . A multifluid scheme is more 
appropriate when explicitly considering chemical reactions and 
the associated change of mass and momentum, as, for example, 
in the dynamics of shock waves in the ISM (Flower, Pineau des 
Forets & Hartquist 1985, Draine 1986) or the collapse of mag- 
netized molecular clouds (Ciolek & Mouschovias 1993). 

In the following, we make the hypothesis that gas is not sub- 
ject to ionization, recombination, or other chemical reactions 
that contribute to the transfer of momentum from one fluid to 
the other In other words, we restrict ourselves to the study of 
nonreacting MHD flows. 

2.5. Equations for the mean fluid 

Neglecting the rate of change of density and momentum of each 
species by chemical reactions (S ss' = rjj' = 0), the equations 
of continuity and momentum for the mean fluid can be easily 
derived summing the equations ( fT9] l and ( |20l i over the species s, 



Dp 

— H-pV-U^O, 
Dt ^ 



(25) 
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DV 



-VP- 



pW + - X B, 

c 



(26) 



where D/Dt = d/dt + U ■ V is the convective derivative for the 
mean fluid. This set of equations must be coupled to Maxwell's 
equations 

'c~dt' 



VxE 



(27) 



VxB = — J+- — , 
c c at 



and 
V B 



0. 



(28) 



(29) 



For a self-gravitating medium, one must also add Poisson's 
equation, 



V^n^ = AnGp. 



2.6. Evolution equation for fhe eiectric current 



(30) 



The mean velocities Uj of each fluid can be expressed in terms 
of U, J, and J^ as 



U 



— (3 -3d), 
q-p 



u+ = U - — J + —3d 

q-p \ Pc 



U-^(^J.-^ 



q-p \p. 



-Pn 



(31) 
(32) 
(33) 



The equation for the electric current J is obtained adding to- 
gether the equations of momentum of each species (Eq.l20b mul- 
tiplied by qs/nis, substituting the expression for the friction force 
(Eq. |2T]i, and eliminating the velocities of each species with the 
help of Eqs. dSTTl-dSSll. The result is: 

^ + V ■ (UJ + JU) = -^V(P_ - eP+) 
ot ni- 

,(iil^(E.HxB 

+±- 

q-Pn 



l-(l+e)^ 
P 

Jrf X B + — 



T „ (1 + e)Q'-+ + a_„ 
J X B J 



P- 



ea+n . 



(34) 
rn^cp Pc 

Equation (l34l i generalizes the evolution equation for the elec- 
tric current derived by, e.g., Rossi & Olbert (1970) and 
Greene (1973), by including explicitly the effects of particle col- 
lisions. 



2.7. Evolution equation for the ambipolar current 

To obtain the equation for the ambipolar current, it is conve- 
nient to subtract the equations for the charged particles from the 
equation for the neutrals in Eq. ( l20b . Expressing the velocities 
of the single species in terms of the global fluid properties with 
Eqs. Olll, (Hill, (EB, ^ and ^, we obtain: 



u> 



H-V- 






^(UJrf + JrfU) 

p 

q-Pn T D "-n T 
- J X B -H J 



m^cp 



P- 



id. 



where we have defined 

_ p„VP-pVP„ 
\j„ = . 



^(UJ + JU) 

p 

(35) 



(36) 



Equation ( |35] ) can be regarded formally as an equation for 
the evolution of the ambipolar current J^, analogous to Eq. (l34b 
for the electric current. The quantity G„ in Eq. (l35T l has the di- 
mensions of a pressure gradient and is a source term for the ion- 
neutral drift that originates from density and/or temperature dif- 
ferences in the two components. 

2.8. Evolution equation for the magnetic field 

Equations ( l34l i. (l35T l. Faraday's law (Eq.lZTb and Ampere's law 
(Eq. l28T l constitute a set of evolution equations for J, 3d, B, and 
E. However, since the time variations of these four quantities 
occur on different time scales, a considerable simplification and 
a deeper physical insight are possible depending on the problem 
at hand. This procedure of separating the timescales has long 
been known, but is often formulated in a qualitative and ad hoc 
manner, as pointed out by Vasyliunas (2005). In this section, we 
examine carefully the timescales associated with each evolution 
equation. 

Let the characteristic length and timescales for the varia- 
tion of fluid quantities be equal to { and t ^ f/U, respectively, 
where U - |U| is the modulus of the average velocity defined by 
Eq. (O. Consider now Eq. (|34| ) for J. The time derivative and the 
convective derivative of J on the left-hand side are both of order 
~ J/t, where 7 = |J|. With the help of Eqs. ^ and ^, it can 
be easily shown that these two terms are much smaller than the 
term proportional to E -i- (U/c) x B on the right-hand side, if 



^»^„i„ = (l+e)-i/2^_, 



where 



A^ = 



m±c 



{W±p±)"^ ^± 

is the so-called inertial length and 



(37) 



(38) 



(39) 



is \hs plasma frequency in the charged particles. It is easily seen 
that ^n,in is equal to the inertial length of the charged particles of 
smallest mass. Notice that the condition for charge neutrality. 



( » Ay),± 



a)± 



(40) 



where /Id,± is the Debye length in the charged particles, is auto- 
matically satisfied if Eq. ( |37l ) is satisfied, because fl± ^ c. When 
the condition dSTJ i is satisfied, the left-hand side of Eq. ( |34] | can 
be neglected, and Eq. ( |34| ) becomes an instantaneous equation 
for E in terms of J, J,;, and B. 

Consider now Eq. ( l35l l. The time derivative and the convec- 
tive derivative of 3d on the left-hand side can be neglected with 
respect to the last term on the right-hand side, proportional to J^, 
if 



T » TcoU = 



(1 + e)p+p„ 



1+e 



(41) 



ia+„ + a-„)p pjAD 
where we have defined the ambipolar diffusion drag coefficient 



Tad = r+« + ey- 



(42) 
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Table 1. Characteristic length and timescales required for the simplification of the evolution equations for the electric and ambipolar 
currents, Eq. d34] i and Eq. ( [35l l. Here the neutral density n„ (in cm"^ is the density of atomic hydrogen and fi is the mean molecular 
weight in units of the hydrogen mass. The electron and ion densities n^ and «, are in cm"^, the magnetic field B is in G, the grain 
radius rg is in units of /im, the interior density of grains is assumed to be equal to 2 g cm"-'. 



composition 


T 


f . 


Tcoll 


Toy,- 




(K) 


(cm) 


(s) 


(s) 


H2, HCO+, e 


10 


5.31 X IQ'nJ" 


1.79xlO'"(M««r' 


5.69 X 10-*'fi-' 


H2, H+, e 


10 


" 


2.51 X lO'Ounn)-' 


" 


H2, H+, e 


10 


" 


2.56 X lO'OuHn)-' 


" 


H, H+, e 


10* 


" 


1.01 X lO'^OunJ-' 


" 


Hz, HCO+, g- 


10 


1.23 X 10**n"i'^ 


2.83 X 10'VjOun„r' 


5.23 X lOV^fi-' 


^2,8\8- 


10 


3.60xl0iVfVi/2 


" 


" 



Thus, the characteristic timescale of the system must be larger 
than an effective collision timescale. Similarly, the time and con- 
vective derivatives of (p„/p)J are negligible with respect to the 
term proportional to J x B on the right-hand side if 



T » T„vr = li^-l 



where 



Q± 



q±B 

m±c 



(43) 



(44) 



is the cyclotron frequency of charged particles (negative for neg- 
atively charged particles) and B - |B|. If conditions (1411 1 and ( l43T l 
are satisfied, the terms containing the space and time derivatives 
of Jd and J can be neglected in Eq. ( l35l l. 

Finally, a comparison of Faraday's law (Eq. IZTT i and 
Ampere's law (Eq.|28Tl, shows that the time derivative of E (the 
displacement current) in Eq. (l28T l is of order {I I erf- x (U/c)^ 
with respect to the curl of B. Thus, in the following, we will 
neglect the displacement current in all applications of Ampere's 
law to ISM conditions. 

Summarizing, the different levels of approximations adopted 
are as follows: 

1. The time and convective derivatives of J are negligible 
on length scales larger than ^^in, the inertial length of 
charged particles of the smallest mass (the so-called large- 
scale plasma limit), regardless of the frequency of collisions 
among particles (see Vasyliunas 2005). 

2. The time and convective derivatives of J^ are negligible on 
timescales larger than TcoU, the typical timescale of collisions 
between charged and neutral particles, and Tgyr, the inverse 
of the cyclotron frequency |Q-| ' of negative charges (the 
collision-dominated plasma limit, see e.g. Dungey 1958). 

3. The time derivative of E (the displacement current) is neg- 
ligible if the characteristic speed of the system satisfies 
(U/c)^ « 1 (nonrelativistic or weakly- relativistic regime, 
see e.g.. Cowling 1957). 

In Table [1] we list the numerical values of these length and 
timescales for several chemical compositions and for physi- 
cal conditions typical of molecular clouds and protostellar jets, 
using the values of the collisional coefficients determined in 
Paper II. 

When the conditions expressed by Eqs. d37T i, Wii and (|43T l 
are satisfied, Eq. (|35] | becomes an instantaneous expression for 
the ambipolar current J^, 



1-Pc 



(a+„ + a-n)m- 



m-a-„ p„ 

J + —(J X B) - G„ 

1-P- cp 



(45) 



which, substituted in Eq. ( [34l i, neglecting the time and convec- 
tive derivative of J neglected, gives the instantaneous expression 
forE, 



„ U „ m_ 

E-H-xB= 

c (1 + e)q-p- 



\ a+„ + a-„ 



Pn 



(a+„ + a^„)pc 



-G„xB-N ^J"«^"« 



cp 



a.„ + a^,, 



q-cp-p 

2 



q-p- 



1 - f [a+n- a-„ 
P- + 7 IP" 



JxB 



a+„ + a^n 



(46) 



relating the electric field E to J and B (one of the forms of the 
so-called generalized Ohm's equation). 

Inserting Eq. ( |46] l in Faraday's law and using Ampere's law 
without displacement current, we obtain the complete evolution 
equation for B in the reference system of the average fluid 

^ - V x (U x B) = — V(P_ - eP+) X Vp_ + V X (^GJ 
at p- 

-V X (^G„ X b) - V X {^B X [(V X B) X B] 



-Vx 
where 



(V X B) X B 



V X (770V X B), 



AfB = 



m^c 



(1 + e)q+p+' 
ea+„ - a-, 



Afg =/^B 



Pn 

m^ — 



a+„ + a- 
P 



p \a+„ + a-„ 
B^ 



'7AD 



^H 



pI 



P+PYAD 

p„B^ 



Anp^ \ a+„ + a^n j 4-7:p^p+yAD ' 



m+cB 
4nq+p 



(1-.)+!^^" 



1 / m+c 
?7o = — I I I a-+ + 






yo / m+c 
a+„ + a-n I 4-n\ q+ 



An \q+p 
and we have defined the Ohmic drag coefficient 



^n y+nj-n 
^+ \ TAD 



(47) 
(48) 

(49) 
(50) 

(51) 
(52) 

(53) 
(54) 
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Table 2. Momentum transfer rate coefficients (crv)sj' and coUi- 
sional drag coeflicients yss' for various chemical compositions 
(see Paper II for analytical formulae and numerical results). For 
collisions between charged particles, the value of the Coulomb 
logarithm has been taken equal to 17.9 and 20.6 for molecular 
clouds and jets, respectively. 



species 


ss' 




<crv>i.,' 


(cm3 S-' g-') 




molecular clouds (T - 


: 10 K) 


HCO+, 


H2 




1.73 X 10-' 


3.33 X 10" 


H+,H2 






1.99 X 10-' 


2.38 X 10'* 


H+,H2 






1.17 X 10-' 


2.33 X lO'-* 


H+,He 






1.43 X 10-' 


1.71 X 10'-* 


e,H2 






1.97 X 10-' 


5.88 X 10"* 


e,He 






1.35 X 10-' 


2.02 X 10"* 


g" (r, = 


: 10 A), 


H2 


1.77 X 10-' 


2.11 X 10" 


8^ ('> = 


: 1 yUm), 


H2 


1.77 X 10-3 


2.11 X 10** 


S" (r, = 


: 10 A), 


He 


1.25 X 10-' 


2.11 X 10" 


g" (rs = 


: 1 yUm), 


He 


1.25 X 10-3 


2.11 X lO'* 


e,H+ 






2.06 


1.23 X lO^-* 




protostellar jets (T = 


10* K) 


H+,H 






2.33 X 10-'* 


5.96 X 10'^ 


H+, He 






8.97 X 10-'" 


1.07 X 10'* 


e,H 






1.44 X 10-^ 


8.61 X 10'* 


s,H+ 






7.48 X 10-5 


4.47 X 10" 



In the limit of very high and very low ionization, the ambipo- 
lar. Hall and Ohm resistivities reduce to the expressions given in 
Appendix A. Equation (l47l i is the fundamental evolution equa- 
tion for B of the problem under the given assumptions. Let us 
now analyze in detail each term of Eq. Wh and evaluate the as- 
sociated diffusion timescale. 



2.9. Ambipolar diffusion 

We consider first the fourth term on the right-hand side of 
Eq. (|47T l, representing the evolution of magnetic field by ion- 
neutral drift or ambipolar diffusion (Mestel & Spitzer 1956). 
The associated resistivity coefficient t/ad is inversely propor- 
tional to the drag coefficient jad, defined by Eq. ( l42l i. Table |2] 
shows that the contribution of electron-neutrals collisions to jad 
is much smaller than the contribution of ion-neutrals collisions: 
for e <K 1, 7ad ~ J+n', conversely, if the negative charge is car- 
ried by dust grains (e » 1), their contribution to ^ad is dominant 
(see e.g., Nakano et al. 2002), and jad ~ ey-«- The multifluid 
expression for the resistivity (see Appendix B) leads us to un- 
derestimate the value of 77 ad by a factor ^j^, and is therefore valid 
only for weakly- ionized gas. In a protostellar jet, for example, 
where the neutral ion fraction can be as low as ^„ =s 0.5 (see 
Sect. |6]l, the actual value of //ad is a factor of ~ 4 and higher 
than that obtained with a multifluid approach. From Eq. ( l47b . 
the timescale of ambipolar diffusion is 



fAD ~ — I Pr AD, 

^n \vAi 

where va - BliA-npYI'^ 
fluid. 



(55) 



is the Alfven velocity in the average 



2.10. Biermann's "battery" 

The first two terms on the right-hand side of Eq. (l47T l rep- 
resent the generation of magnetic fields by electric currents 



produced by pressure gradients in the charged fluid, the so- 
called Biermann's "battery" (Biermann 1950, Schliiter & 
Biermann 1950). The first term represents a pure plasma process, 
whereas the second involves collisions with neutrals. Generation 
of seed magnetic fields by Biermann's battery processes are pos- 
sible only if the density and temperature gradients of the fluid 
components are not aligned, a condition that may occur in stel- 
lar interiors but is unlikely in the ISM. It is easy to verify that 
the first term is larger than the second if the negative charge is 
carried by electrons, whereas the opposite is true in the case of 
negatively- charged grains. 

2.11. Diamagnetic current 

The third term represents the effect of a diamagnetic current, 
proportional to G„ x B, and perpendicular to the magnetic field. 
In the presence of a magnetic field and a pressure gradient, the 
charges move across regions of different densities producing a 
drift current perpendicular to the field and the pressure gradi- 
ent, that tends to reduce the strength of the field in the plasma. 
Diamagnetic currents are generally negligible in the ISM except 
perhaps near cloud's boundariefl There are, however, situations 
in which diamagnetic currents may contribute significantly to 
the dissipation of the magnetic field and the heating of the gas. 
If the temperatures T, and the mass fractions ^j of the fluid com- 
ponents are spatially constant, one finds 



G„ = ^„(fl2 - fl2)Vp, 



(56) 



showing that G„ vanishes both in the limit of very high (^„ — > 0) 
and very small ionization {a„ — > a). Figure [1] shows the value 
of G„ calculated for some gas mixtures considered in Paper II, 
assuming uniform temperature and composition. It is clear from 
this figure that for intermediate values of the ion fraction (^+ ^ 
0.2-0.8, depending on the ISM) composition, the magnitude of 
G„ is ~ 40% of the pressure gradient, and, therefore, is not neg- 
ligible a priori. With G„ given by Eq. (|56l l, the diamagnetic dif- 
fusion timescale reads 



.D.^ 



e^ 



1 2 /''''AD 



& w 



independent on the intensity of the magnetic field. 



(57) 



212. Hall diffusion 

The three-fluid scheme reveals that the Hall resistivity rjn, de- 
fined by Eq. (|52] |. is made of two terms. The first term van- 
ishes when the charged particles have equal mass-to-charge ra- 
tios {e - 1), whereas the second vanishes when the particles have 
equal collisional coefficients (a^„ - a+„). Only the latter term 
appears in the multifluid approach usually adopted for a weakly- 
ionized gas (see Appendix B). If a-„ + a+„, the multifluid ap- 
proach is valid when the first term inside the square parenthesis 
in the expression of rju is smaller than the second term, of the 
order of ^„/^+, i.e., when the condition 



1 



& \l-e\ 



(58) 



is satisfied. For a plasma of neutrals, ions and electrons (e <K 1), 
the condition ( |58] | is equivalent to ^+/^„ <k 1, largely satisfied 



' In the geomagnetic context, the Chapman-Ferraro current flowing 
at the interface between the solar wind and the Earth's magnetopause is 
a diamagnetic current (see e.g., Parks 1991). 
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and is independent of the density and the intensity of the mag 
netic field. 



3. Advection-diffusion equation for the magnetic 
flux 

A more compact form of Eq. ( |47] i can be obtained separating the 
poloidal and toroidal components of the velocity and the mag- 
netic field. 



UsU/,+U^, 



B = Bp -H B^. 



(64) 



In the following, for simplicity, we will consider only axially- 
symmetric systems {djdip - 0). In a system of cylindrical coor- 
dinated (R, z, if), we introduce a vector potential 



^(R,z) 



0.10 



1.00 



HR,z)= ^^^ e. 


' 


such that 


B;, = V X A = V X 





Fig. 1. Function G„ normalized to the pressure gradient in the 
average fluid VP as function of the mass fraction of positive 
charges ^+. Dotted curve: electrons, HCO^, H2; short-dashed: 
electrons, H^, H2. long-dashed curve: electrons, H^, H2; dot- ^{R,z) - 2nRB^{R,z}. 



1 
and a scalar function ^ defined as 



50), 



5<D, 



-Tz'^ ^ m'' 



(65) 



(66) 



(67) 



dashed curve: electrons, H+, H. 

in molecular clouds and, marginally, in protostellar jets. In this 
case, Eq. ( |52l ) reduces to the same expression obtained with the 
multifluid approach (see Appendix B) 



m~ j- 



^n i m+cB 



Anq+p 



(59) 



Equation ( |47l ) can be separated in one equation for B^, and 
one equation for B^^. Expressing B^ as function of cD according 
to the definition (|66] |. the induction equation for B^ can be "un- 
curled" and expressed in scalar form as 

^+U,.VcI.= ^G„.VO 
dt '^ P 



riAD 



independent on the coUisional coefficients. The associated Hall 
diffusion timescale is 



(InRB)^ 

m 



[.S(<1))V(1) -I- ^PV*?] ■ VO 



fe 



^n \vai 



InRB 



J-(*I',<D) + 77o>S(<D), 



(68) 



(60) 



where B - |B|, S and J' are Stokes and Jacobi operators defined 

as 



If the negative charge is carried by ISM grains of radius rv, 
the multifluid expression for the Hall resistivity given by Eq. ( |59] | 
is a good approximation of the correct expression Eq. ( |52] | 
only when the ion fraction is extremely low, ^+/^„ <K e"' « 
lO^^^irg/fixny^. For larger mass fractions of positive charges, 
typical of ISM conditions, the Hall diffusivity becomes indepen- 
dent of the ion fraction. 



5^ 5^_j_50 
^ ^ ^ M^ ^ W RM' 

and 

^^ -^ dR dz dz dR' 



(69) 



(70) 



m 



m^cB 

A-nqp 



(61) 



In this case, the associated Hall difi'usion time scale becomes 

fH * (£j l^-l- (62) 

In this case, the multifluid expression for rju given by Eq. (|59T l 
leads to a serious underestimate of the actual value of the Hall 
resistivity. 

2.13. Ohmic diffusion 

The Ohmic diffiision timescale is 



respectively. 

In the particular case of zero toroidal field, Eq. ( |68] ) reduces 
to the simple expression 



^ + U^ ■ V(D = ^G„ ■ VO + (;7AD + /7o)>S(<D), 
ot P 



(71) 



fo^ - 



va/ epro 



(63) 



In this reduced form, and neglecting the diamagnetic term, 
Eq. ( TTTI ) is at the basis of many studies of magnetic flux evo- 
lution by ambipolar and Ohmic diffusion in molecular clouds 
(see e.g., Desch & Mouschovias 2001). 

Ambipolar diffusion produces a time variation of the mag- 
netic flux enclosed in a given region, at a rate depending on 
the ionization fraction and the strength of the poloidal and 
toroidal field. Consider, for example, the realistic case of a ver- 
tical magnetic field concentrated toward the z axis, for which 
<D = Q){R) oc R" with < n < 2. Since >S((D) < for this con- 
figuration, ambipolar diff'usion will drive the field outward of 
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the region under consideration, thus reducing the magnetic flux 
locally. A toroidal field will speed up (slow down) this process 
of field diff'usion, depending on whether the quantity ^' increases 
(decreases) outward, or in other words, whether the toroidal field 
decreases slower (faster) thanT? '. A toroidal field B^ oc 7? ' has 
no effect on the rate of flux loss driven by ambipolar diffusion. 

The effect of the Ohmic diffusion term is to decrease (in- 
crease) the magnetic flux enclosed in a given region if S(<i>) is 
negative (positive). For example, for a magnetic flux O - 0(/?) oc 
R", the function iS((B) is positive (negative) if n > 2 (n < 2), or, in 
other words, if the magnetic field strength increases (decreases) 
outwards. The presence of a toroidal field has no effect on the 
variation of the magnetic flux by Ohmic diffusion. 



Table 3. Ambipolar diffusion coefficients ^ad and x for colli- 
sions with H2 in molecular clouds. The values of ;t' ^re obtained 
withC = 3.9 X 10-1^ g'/2cm-3/2^ndC = l.SxlQ-" g^'^ cw.-^'^ 
for non-depleted and depleted clouds, respectively. 



composition 


Tad 

(cm3 s-i g-i) 


X 


H2, HCO+, e 


undepleted cloud 

3.33 X 10'-'' 


1.0 


H,, HCO+, g- 


(rj = 10 A) 6.97 x 10" 


2.1 


H,, H+, e 


depleted cloud 

2.38 x 10'" 


2.8 


H,, H+, e 


2.33 X 10'" 


2.7 


H2, H+, g- (r, 


= 10 A) 1.29x10'-^ 


15 



4. Rate of energy generation 

The energy generated by the friction between streaming par- 
ticles is a significant source of heating for molecular clouds 
(Scalo 1977, Lizano & Shu 1987) and protostellar winds and jets 
(Ruden, Glassgold & Shu 1990, Safier 1993). The energy pro- 
duced heats the bulk of the gas, increasing its pressure and allow- 
ing chemical reactions to proceed at a faster rate (Flower, Pineau 
des Forets & Hartquist 1985, Pineau des Forets et al. 1986). 

In the absence of external forces, the rate of energy genera- 
tion (per unit time and unit volume) is the sum of the work done 
by friction forces on the three species. 



W = ^u,-^F,,, 



= a-„|u- - u,/ + a+„\u+ - u„f + o-^+lu- - u+f. 



Substituting in Eq. ( |72] | the expressions for the drift velocities in 
terms of J and J,/ obtained from Eqs. ([3TTl-(l33]l, and eliminating 
Jd with the help of Eq. ( l45T l. we obtain 



W = 



P+PnTAD 



Pi, 

Anp 



[B x (V x B)] -I- G„ 



+ ^|VxB|2 
A-n 



(73) 



(see also Braginskii 1965). The first term on the right-hand side 
represents the dissipation of drift motions between the plasma 
and the neutrals driven by the Lorentz force or pressure gradi- 
ents, whereas the second term represents the heating of the gas 
due to the dissipation of the electric current. The ratio of the sec- 
ond to the first term inside the modulus can be estimated with 
the help of Eq. 



47rp|G„| 



p„|B X (V X B)| 



(74) 



Equation ( |74| ) shows that the contribution of pressure effects to 
the heating of the gas is important only in weakly-magnetized 
and/or hot plasmas with significant ionization. As shown in 
Fig.[T] the pressure-driven drifts, represented by G„ in Eq. ( l73T l 
are significant only at intermediate levels of ionization. 
If one sets G„ - 0, the heating rate is 



W ■ 



1 






r]AD- 



|>S(<D)V<D + ^FV^Fp + \J(^', (D)|2 



|V(&p + iv^pp 



■ |V*Pft] . 



(75) 



In the particular case of zero toroidal field, Eq. dTST l reduces to 
the simple expression 



* = ^^i«*)P- 



(76) 



5. Magnetic field diffusion in molecular clouds 

To illustrate the formalism developed in the previous sections, 
we analyze the problem of determining the diffusion of the mag- 
netic field in two different astrophysical environments: cold, 
weakly-ionized molecular cloud cores, and hot, mildly-ionized 
protostellar jets. In this section, we analyze the case of molecu- 
lar clouds, while protostellar jets are considered in Sect.|6] 

In dense molecular clouds with T ^ 10 K, n„ « 10^- 
10'' cm^^ «/ X 10" '-10^2 cm"^ and B ^ 10"^ G, electrons 
are in general the main carriers of negative charge (Nakano et 
al. 2002). In these conditions, the ion density is determined, in 
a three-fluid scheme, by a balance of cosmic -ray ionization and 
ion recombination, resulting in the simple relation 



(72) p^^Cpl 



1/2 



(77) 



In nondepleted clouds, the dominant ionic species are metal ions 
(like Mg+, Na+, etc.), molecular ions (like HCO^), or H+ and 
H^ ions, the simple Eq. ( |77] ) is in good agreement with the 
results of detailed chemical models (Nakano et al. 2002) with 
C « 3.9 X 10"'^ g'^^ cm"^''^. In completely depleted cores, 
where all metal species are frozen onto grains, the dominant 
ions are H:^, H^ and their deuterated analogues. In this case, 
the ion density can be approximated by Eq. ( TTTI i with C ~ 
1 .5 X 10"'^ g'^^ cm"^''^, although the actual ion density increases 
with increasing density slightly less steeply thanp'^^ (Walmsley 
et al. 2004). In Table|3] we list the values of the ambipolar diffu- 
sion coefficients for collisions with H2 appropriate for different 
chemical compositions and degrees of depletion (we ignore the 
10% correction due to collisions with He, see Paper II). 

It can be easily verified from Table [1] that the conditions 
for the validity of our three-fluid MHD approach, expressed by 
Eqs. (IJTT i. dJU and (l43T l are largely satisfied. However, a compli- 
cation arises at high densities (n„ x 10"-10'^ cm"-') and rela- 
tively small scales characterizing the so-called decoupling stage 
of matter and of magnetic fields, where the electric charge is 
mostly carried by dust grains with ng- » n^+ ^ 10""n„ (Desch 
& Mouschovias 2001, Nakano et al. 2002). A comparison be- 
tween the typical size of a molecular cloud during the dynamical 
coflapse, ^ s; 10 (njcm'^y^^^ pc (Nakano et al. 2002), with the 
value of ^n,in listed in Table [T] for a plasma of negatively- and 
positively-charged grains with radius rg shows that the condition 
^ » ^min is satisfied for ng-/n„ » 10"'^(r^//im)-'. This con- 
dition is violated if grain coagulation during cloud collapse is 
rapid, and the mean grain radius can reach yum-size dimensions 
(Flower, Pineau des Forets & Walmsley 2005). Thus, an accu- 
rate analysis of the decoupling stage must take into account the 
evolution, during the collapse, of their mean radius to ensure that 
conditions ( |37] l is satisfied at all times. Otherwise, the evolution 
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equation for the field, Eq. ( |47T i must be replaced by the evolution 
equation for the current, Eq. ( [34] i. with B given as function of J 
by Faraday's equation (Eq.lZTJ). 

The most important diffusive process in molecular clouds is 
ambipolar diffusion, occurring on a timescale 



f AD =» — TADP, 



(78) 



(Mestel & Spitzer 1956). It is easy to verify that the timescales 
of the Biermann's "battery" and diamagnetic diffusion are too 
long in molecular cloud conditions. In typical molecular cloud 
conditions, if the negative charge is carried by electrons, the 
timescale for this process is of the order of ~ 10'*(^„/^,) times 
the ambipolar diffusion timescale, and therefore can be safely 
ignored. Also, Biermann's battery terms identically vanish when 
the charge is carried by positively- and negatively-charged grains 
of equal properties. However, the Biermann's battery may lead to 
the generation of seed magnetic fields at the boundary between 
a dense atomic or molecular cloud and a hot diffuse external 
medium where gradients in the density, temperature or chemi- 
cal composition can be strong (Lazarian 1992). For diamagnetic 
effects, a comparison of Eqs. ( |55] l and ( |57l ) shows that 



Jr. 

fAD 



a^ - at 



(79) 



a large quantity in weakly-ionized molecular clouds where (a^ - 
al)/a^ ss 2^i{m„lmi) <K 1 and va > a. 

In general, however, the Hall term in Eq.(|68]l cannot be ne- 
glected with respect to ambipolar diffusion (Wardle 1999). If 
the low-ionization condition ( fSSl l is satisfied, we obtain from 
Eqs. (|55]l and ^ 



tAD P7ad ' 



(80) 



independent on the ion fraction ( Wardle & Ng 1999). The quan- 
tity on the right-hand side of Eq. dSOl l is called the ion's Hall 
parameter (see Appendix B). Inserting the values of jab listed 
in Table [3] we obtain 



fAD \yuG/ \cm ^ I 



(81) 



for depleted and undepleted clouds alike. Thus, for densities 
n„ ~ 10*' (BIjiGY^ cmT^ or larger, the influence of Hall diffusion 
competes with ambipolar diffusion in determining the evolution 
of the magnetic field in a molecular cloud. The importance of 
Hall diffusion is even larger when the negative charge is carried 
by dust grains and the condition (|58] | is not satisfied. In this case, 
the ratio of Hall and ambipolar diffusion timescales is 

iiL . k (B-\ , iQ-^^k (I.\ (JhLS' (^r , (82) 

fAD ^i \py AD I ^i \/^G/\cm"^/ \^imj 

and may become of order unity or less at sufficiently high densi- 
ties. The effects of Hall diffusion have yet to be incorporated 
in realistic models of cloud collapse. The importance of the 
Hall diffusion for molecular clouds (and circumstellar disks) has 
been stressed by Wardle and coworkers (Wardle 1999, Wardle & 
Ng 1999, Pandey & Wardle 2007). The evolution of the cloud's 
magnetic field under Hall diffusion is a complex nonlinear pro- 
cess, strongly coupling the poloidal and toroidal components 
to each other (Goldreich & Reisenegger 1992, Wardle 1999), 



which has yet to be incorporated in realistic models of cloud 
collapse. 

Finally, the Ohmic resistivity is often assumed to be 
dominated by electron-neutral collisions (e.g., Desch & 
Mouschovias 2001), because in the multifluid treatment, where 
only collisions of charged particles with the neutrals are con- 
sidered, electrons have the largest Hall parameter /?(.„ (see 
Appendix B). This is hardly the case in molecular cloud condi- 
tions, where the large value of the Coulomb cross section com- 
pared to the polarization cross section (see Paper II) makes colli- 
sions between charged particles nonnegligible or even dominant 
with respect to collisions with neutrals if the ionization fraction 
is sufficiently high. To illustrate this point, we show in Fig.|2]the 
value of the critical ion fraction above which the Ohmic resistiv- 
ity is dominated by collisions between charged particles, corre- 
sponding to the condition q'_+ > a+„Q'_n/(Q'+„+Q'_„)(seeEq.|53]l. 
The critical value of the mass fraction of positive charges ^+/^„ 
is shown as a function of the neutral density, assuming a tem- 
perature appropriate for a collapsing molecular cloud core fol- 
lowing Tohline (1982), for various chemical compositions of 
the gas, with the collisional rate coefficients given in Paper II. 
For typical conditions of dense cores, n„ ^ 10"^-10^ cm"-', 
electron-ion collisions are dominant over electron-neutral colli- 
sions for most chemical compositions, and the multifluid expres- 
sion of 770 underestimates the value of the Ohmic resistivity (see 
Appendix B). However, this limitation of the multifluid approach 
has little consequence for the evolution of molecular clouds and 
cores because ambipolar diffusion overwhelms Ohmic dissipa- 
tion by several orders of magnitude. Interestingly, for the condi- 
tions of the so-called decoupling stage, when the magnetic field 
is expected to decouple from the matter in a collapsing cloud, 
with n„ X 10"-10'^ cm"', the charge is carried mostly by dust 
grains with pg- ^ pg+ (Desch & Mouschovias 2001, Nakano et 
al. 2002). Under these conditions, depending on the grain radius, 
and the fraction of charged grains, the contribution of grain-grain 
collisions to Ohmic dissipation may even dominate over grain- 
neutral collisions, a possibility that has never been considered in 
calculations of the decoupling stage of star formation. 

5.1. A toy-model molecular cloud 

To estimate the timescale of magnetic flux redistribution and the 
associated heating rate in the typical conditions of a molecular 
cloud core, we consider a toy model for an isothermal, cylin- 
drical gas cloud in equilibrium under the effect of gravitational, 
magnetic and pressure forces. The hydrostatic equilibrium of the 
cloud is described by the Eq. (|26] | with d/dt = and U = 0, cou- 
pled to the Poisson equation, Eq. ( l30b . A simple solution satis- 
fying these equations was found by Nakamura et al. (1993), and 
reads 

p(x) = 77-^, (83) 



(l+x2)2' 



B« = 0, B.^Bq 



1 +(1-^0^2 



,21 1/2 



(1 H-x2)3 



(x 

^f ~ ^0— JTTJJ 



,(84) 



where po and Bq are the central values of density and magnetic 
field, ^ is a parameter with value < ^ < 1 measuring the rel- 
ative importance of the toroidal and longitudinal components of 
the magnetic field, and x - R/Rq is the ratio of the radial coor- 
dinate R and the characteristic radial scale Rq. The parameter ^ 
is related to Rq, a, po and Bq by the relation 

^^'nGpoRl 



-e-4 



2-C 



\VAfi 



2fl2 



- 1 



(85) 
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Fig. 2. Importance of collisions between charged particles in 
molecular clouds. In the area above the curves, collisions 
between charged particles dominate over collisions between 
charged particles and neutrals in the expression of the Ohmic 
diffusivity //o- The chemical composition of the gas (neutrals, 
positive, and negative species) is: H2, H^, e {dotted curve); H2, 
H3, e (short-dashed curve); H2, HCO^, e (long-dashed curve); 
H2, HCO^, g^ (dot-short dashed curve); H2, g^, g^ (dot-long 
dashed curves, with rg = 50 A and r^ = 100 A). The tempera- 
ture is shown by the solid curve. 



where G is the gravitational constant and v^.o = Bo/(47rpo)^^^ is 
the Alfven speed on the axis of the cloud. 

If the ionization fraction can be approximated by Eq. ( iTTb . 
the time scale for ambipolar diffusion in a magnetically- 
subcritical cloud from Eq. ( |68] l is 



tgfl \2a) 



R^^\^(nGpoRl 



2fl2 



-1 



Tc\(x, 



where 



X ^ 



7adC 
(87rG)i/2 



(86) 



(87) 



is nondimensional parameter of order unity or larger (Shu 1983, 
Gain & Shu 1993a,b), fff^o = (InGpo)-^'^ is the free-fall 
timescale on the axis of the cloud, and Tci(x, () is the function 
shown in Fig. [3] In Eq. ( l86b . the approximate equality is valid 
for strongly-magnetized clouds where v^fi » a. 

For this particular model, it is easy to check that the ambipo- 
lar diffusion term on the right-hand side of Eq. ( |68] l is negative 
for all values of ^, and therefore the effect of ambipolar diffu- 
sion always decreases the flux contained in fluxtubes close to the 
cloud's axis. The time scale of the process depends, however, on 
the strength of the toroidal field. As shown by Fig. [3] the diffu- 
sion time of the magnetic flux increases with distance from the 
cloud's axis and with increasing strength of the toroidal compo- 
nent of the field, up to a factor of ~ 2 on the cloud's axis. 

The relevant collisional rate coefficients for molecular cloud 
conditions computed from the expressions derived in Paper II 
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Fig. 3. Timescale for ambipolar diffusion Tci in the model cloud 
described in Sect. (15.1b as a function of the distance from the 
axis X and the relative strength of the toroidal magnetic field. 
The solid curve is for ^ = 0, the dotted curve for ^ = 1 . 

are listed in Table |2] Figure [3] shows that the evolution of the 
magnetic flux (and therefore of the whole cloud) driven by am- 
bipolar diffusion occurs faster in the central regions of the cloud, 
where Tad is larger than the local free-fall time by a factor x- If 
the innermost parts of the cloud are strongly depleted, as often 
observed, (see e.g., Caselli et al. 1999), the values of the diffu- 
sion coefficient yAD and the parameter;^ listed in Table |3] show 
that the evolution of the magnetic fields occurs on a timescale 
of ~ 3fff_o if the dominant ions are H^ or H^ and the negative 
charge is carried by electrons, or ~ 15fff_o or larger if the neg- 
ative charge is carried by a large number of very small grains. 
An accurate knowledge of the chemical composition of the cen- 
tral regions of a molecular cloud core is thus necessary before 
conducting any evaluation of the core's evolution timescale, if 
ambipolar diffusion is the driving agent. 
The ambipolar diffusion heating is 



WAn = 



2po (2a^f(nGpoRl 



X 



2fl2 



Po 



2j^ffF,0 \fff,0 



^0/ I 

: — Wc\(x), 



- 1 



fff.o w^x(x). 



(88) 



where Wc\(x) - x^/(l + x^)^, and, again, the approximate equal- 
ity is valid for strongly magnetized clouds where v^o » a. 
The maximum heating rate is reached at x - 1 / V2, where 
Wci - 4/27. The Ohmic heating is always negligible compared 
to ambipolar diffusion heating in typical molecular cloud condi- 
tions. 

6. Magnetic field diffusion in protostellar jets 

The diffusion of magnetic fields plays an important role in 
jets from young stellar objects, caractherized by typical length 
Lj&t ~ 1 pc or more; typical radius /?jet ~ 10^ AU; total den- 
sity between lO-' and lO'* cm"-'; and temperature T x 10'* K 
(Reipurth & Bally 2001). Recent studies have shown that the 
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ion number fraction «,/(«, + n„) can vary between 0.01 and 0.5 
(Bacciotti & Eisloffel 1999; Bacciotti, EisloflFel & Ray 1999; 
Podio et al. 2006; Hartigan et al. 2007). The bulk velocity in 
protostellar jets is ~ 200 km s"', so that the flow is highly su- 
personic (Mach number ~ 10). Little is known about the magni- 
tude of the magnetic field, which is difficult to determine (Mundt 
et al. 1990; Hartigan et al. 1994; Ray et al. 1997; Hartigan et 
al. 2007). The intensity of the magnetic field in jets is sometimes 
estimated to be in the range 10^-10^ fiG from equipartition con- 
siderations. 

In Table |2] we list the relevant collisional rate coefficients 
for typical conditions of protostellar jets computed in Paper 11. 
For a jet made of H, H^ and e at T = lO'* K, we obtain 
from Table |2] the value of the ambipolar diffusion coefficient 
Tad - 7.01 x 10'^ g"' cm-' s, with e-H collision contribut- 
ing less than 1% to the total. For the same composition and 
temperature, yo is largely dominated by e-H^ collisions, yo - 
4.47 X 10'^ g"' cm-* s. It can be easily verified from Table[T]that 
the conditions for the validity of our three-fluid MHD approach, 
expressed by Eqs. (|37] ). ( |4TI ) and ( |43] ) are largely satisfied. 

6.1. A toy-model protostellar jet 

We apply now the equations derived in Sect. [3] to a toy model 
of a magnetized jet, represented as a cylindrically-symmetric, 
supersonic flow in equilibrium under the effect of centrifugal, 
magnetic, and pressure forces. We consider an infinite cylinder 
with an axis in the z direction, and characterized by a radial scale 
Ro, not necessarily coincident with the radius of the optical jet 
Rj. All the variables are assumed to be independent of z and (p 
(cylindrical and axial symmetry), and the average flow velocity 
U is everywhere parallel to the magnetic field B. We also assume 
Ur - Br = Q and an isothermal equation of state, P - a^p. 

In steady state, and neglecting the gravitational potential 'V, 
Eq. ( |26] l becomes 



U'^ dP B'^ Id 



A particular solution of eq. ( [89] l is given by 



p(x) = 



Po 

1-Hx2' 



with Po = 



Bl 



4;rt/2 ' 

jet 



U^(x) = t/jet, U^(X) = UjetOx, 



B,(x) 



Bo 



(1 -Hji:2)i/2' 



B^(x) = Bo 



ex 



(l-Hx2)l/2' 



(89) 



(90) 



(91) 



(92) 



where x - R/Rq, Aljet = f/jet/a » 1 is the Mach number of the 
flow along the jet, and = {I + 2/M^y^^. Since the observed 

jel 

Mach number in stellar jets is Aljet ~ 10, we have ^ 1. In this 
particular solution, the toroidal component of the velocity (and 
of the magnetic field) increases linearly with R, whereas the lon- 
gitudinal velocity is independent of the distance from the axis. 
An upper limit on the observed rotational velocities U^ in proto- 
stellar jets is ~ 10% of the axial velocity f/jet at axial distances of 
7?jet ^ 102 AU (Bacciotti et al. 2002; Woitas et al. 2005; Coff^ey 
et al. 2007), implying that the solution above cannot be adopted 
to describe an optical jet beyond .^c a; 0.1. If x » 0.1 corresponds 
to the radius of the optical jet /?jet, we must set Rq ^ 107?jet. The 
density and intensity of the longitudinal magnetic field are then 
almost constant across the optically-visible jet, with the inten- 
sity of the toroidal field increasing roughly linearly with distance 
from the jet's axis. 



As in the case of molecular clouds, flie dominant diffusive 
process for the magnetic field in a protostellar jet is ambipolar 
diffusion. Diamagnetic diffusion may also play a significant role 
at the boundary between the jet and the ambient medium. Near 
the jet's axis, however, the timescale of diamagnetic diffusion, 
from Sect. [3] is 



Id 



Bl 



47r(a2 - al)po 



fAD 



M2 
' 'jet 



fAD, 



(93) 



more than two orders of magnitude longer than the ambipolar 
diffusion timescale ^ad- if Aljet ~ 10. Similarly, the timescale 
for Hall diffusion is 



tu 



P7AD 



■ tAD, 



(94) 



In the range of parameters considered, and inserting the values 
of Tad listed in Table|2l we obtain 



iiL.io^(4)(^r'. 

fAD \juG/ \cm ^ / 

The Ohmic dissipation timescale is 



to 



TmA 
ro / 



^„/3i„ tAD, 



(95) 



(96) 



where y6,„ is the Hall parameter of ions (see Appendix B). For the 
assumed values of jet parameters, y6,„ x 10^-10^, and therefore 
tu is much longer than fAD, despite the fact that yo is about four 
orders of magnitude larger than yAo (see Table|2]l. 

For X <s; 1 and Aljet s> 1, the timescale for magnetic flux 
diffusion by ambipolar diffusion is 



l(^i\ 7ADP0Rl 



(97) 



Inserting order of magnitude values for the physical parameters 
in Eq. ( l97b . we obtain, on the jet's axis. 



fAD^6xl04(|^ 



Tad 



Ui 



jet 



102 km s- 



^„/\10'5g-'cm^^ 

-2 



^sjllO^cm-VllO^AuJ 



comparable to the jet's dynamical timescale, 

, -1 



tdyn ~ 



-jet 



f/i 



jet 



104 (^ 

pc 



U^ 



jet 



102 km s 



yr. 



(98) 



(99) 



This important point has already been stressed by Frank et 
al. (1999). However, our expression of t^o (Eq.l98]l differs from 
theirs for two reasons: first, Tad is proportional to ^,/^„ not 
to ^j^n, as in the multifluid description adopted by Frank et 
al. (1999) (see Appendix B); second, their adopted value for the 
collisional rate coefficient (crv),„ is about one order of magni- 
tude smaller than the value we compute in Paper 11; third, the 
curvature radius of the magnetic field R^ may be larger than the 
actual radius of the optical jet 7?jet, as shown by our equilibrium 
model. 

Notice that in this model d(S>/dt near the jet axis is positive 
for any value of Aljet: the strong "hoop stress" due to the toroidal 
fields compresses the plasma and the poloidal magnetic field to- 
ward the jet's axis, where B. increases with time over a timescale 
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/ad- Equation ( |97] l shows that the timescale of ambipolar diffu- 
sion is directly proportional to ion fraction and the rate of mass 
loss in the jet, and inversely proportional to the square of the jet's 
speed. If these quantities are not strong functions of the distance 
from the central star, the results of our toy-model are robust, and 
indicate that the effects of ambipolar diffusion are more impor- 
tant in the outer parts of the jets, where the dynamical timescale 
is longer. In fact, Frank et al. (1999) argued that field diffusion 
already becomes effective at distances as small as a few tens of 
a pc from the central source. 

The ambipolar diffusion heating rate results 



1Vad = 



jet 



M 



Wjeiix, 0), 



(100) 



where Wjet(jc) ~ x if 6* ~ 1. Thus, at any fixed distance from the 
central source, ambipolar diffusion heating is higher in the outer 
parts of the jet rather than close to the jet's axis. The heating rate 
(in physical units) for our toy-model jet as a function of the dis- 
tance from the jet's axis is shown in Fig. H) compared with the 
mechanical heating required to reproduce the observed emission 
of jets, estimated following Shang et al. (2004). It is clear from 
this figure that in weakly-ionized jets, with «,/(«,- -i- n„) « 0.01, 
ambipolar diffusion heating may provide a good fraction of the 
required energy input in the outer layers of the jet. These results 
are in qualitative agreement with those of Garcia et al. (2001), 
who considered in detail the effects of ambipolar diffusion heat- 
ing on the thermal structure of protostellar jets. This work, how- 
ever, was assuming externally a steady-state MHD model for the 
jet structure, not including ambipolar diffusion. The evolution of 
the magnetic field internal to the jet in presence of ambipolar 
diffusion, however, has never been considered self-consistently 
in a MHD model of the jet structure and kinematics. Our work 
indicates instead that this point should be carefully examined in 
the future. 

7. Conclusions 

We have derived a general and self-consistent form of the MHD 
equations for a three-fluid system containing particles of nega- 
tive and positive charge, and neutrals. We have made no assump- 
tion about the ionization degree, or on the particle's masses, and 
we have considered all possible coUisional processes between 
particles of the three fluids. We have shown that the coupled 
equations for the electric current and the drift current are not 
evolution equations above the large-scale plasma limit and the 
collisionally-dominated plasma limit, respectively, and we have 
derived the evolution equation for the magnetic field valid in this 
regime in the reference frame of the average fluid. The result- 
ing expressions of the various resistivity coefficients appearing 
in this equation differ from those usually adopted in star for- 
mation studies by factors nonnegligible when the ion fraction 
of the gas is significant and/or when collisions between charged 
species cannot be ignored. For axially-symmetric systems, we 
have reduced the equation for the evolution of the poloidal field 
to an advection-diffusion equation for the magnetic flux that in- 
cludes the effects of the toroidal field component, and evaluated 
the typical timescales associated with each diffusive process. We 
have also derived accurate expressions for the heating rate pro- 
duced by the collisional dissipation of drift velocities between 
particles of different species driven by magnetic forces or non- 
uniform pressure gradients. 

We have applied our results, combined with the accurate val- 
ues of the collisional rate coefficients calculated in Paper II, to 
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Fig. 4. Ambipolar diffusion heating rate Wad as function of the 
distance from the jet's axis x - R/Rq compared with the mechan- 
ical heating rate empirically determined by Shang et al. (2004) 
for three values of the distance from the central star (top to bot- 
tom: 0.01, 0.1 and 1 pc, dot-dashed lines). The dotted and solid 
curves are for an ion number fraction of «,/(«,■ + n„) - 0.5 and 
0.01, respectively. 



the study of magnetic flux dissipation in molecular clouds and 
protostellar jets. In the former case, our main conclusions are: 

1 . The timescale of ambipolar diffusion in the central regions of 
magnetized clouds is of the order of the free-fall time (~ 1 to 
3 times depending on the degree of depletion) or about one 
order of magnitude longer if the negative charge is carried 
by a population of very small grains. 

2. The presence of a toroidal component affects the timescale of 
ambipolar diffusion by a few factors, showing that a detailed 
knowledge of the magnetic field strength and morphology, 
in addition to the chemical composition, is necessary for an 
accurate estimate of the time scale of magnetic flux loss in 
molecular clouds. 

3. Collisions between positively- and negatively-charged dust 
grains, neglected in previous studies of cloud collapse, may 
significantly increase the value of the Ohmic resistivity at 
the high densities where the magnetic field is expected to 
decouple from the gas, enhancing the rate of the process. 

4. The Hall resistivity can take larger values than previously as- 
sumed, especially when the negative charge is mostly carried 
by dust grains. 

For typical conditions of protostellar jets, characterized by 
higher temperatures and ion fractions, our study shows that: 

1 . The ambipolar diffusion timescale is of the same order of 
magnitude of the dynamical timescale of the jet, in agree- 
ment with previous estimates. Thus, diffusive effects cannot 
be neglected in the study of the MHD properties of proto- 
stellar jets. 

2. A toy-model of a cylindrical jet shows that the hoop stresses 
produced by the toroidal field are likely to force the poloidal 
field to diffuse toward the jet's axis. The resulting evolu- 
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tion of the jet's structure has never been considered in MHD 
models of protoste liar jets. 
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Appendix A: Limiting cases 

A.1. High ionization 



In the limit of high ionization, p„ « 0. The ambipolar diffusion resistivity, be- 
ing proportional to the neutral density, vanishes (;;ad ~ 0), the Hall resistivity 
becoines independent on collisional rate coefficients, 



;7H ~ (1 - e) 



m+cB 



Anq+p 
and the Ohmic resistivity becomes a function of the temperature only, 

^2 - A 



rio 



1 Im+c 



Z+m-m+c 



T<o-v)-+- 



An \ ep+ / 4n-Z_(m_ + m+)Ze^ 

In the high-ionization liinit, the rate of energy generation becomes 

W -- 



^IVxBll 
4;n- 



(A.l) 



(A.2) 



(A.3) 



A.2. Low ionization 

In this limit, p+,p- ~ 0, p„ 
can be simplified as 



: p. The ambipolar diifusion and Hall resistivities 



1 



'?AD~ — 



B^ 



m- 



m+cB 



Anq+p+ 



(1-6)+ 



^n i a+n - a-r, \ 



(A.4) 



(A.5) 



The expression for the Ohmic resistivity cannot be simplified assuming that 
Qr_+, being proportional to the square of the ionization fraction, is negligible with 
respect to a-,, and Q'+„ because of the large value of the Coulomb cross section 
compared to the typical ion-neutral scattering cross section (see Paper II). 

In the low-ionization limit, the rate of energy generation becomes 

W=^|VxBp + ^|Bx(VxB)|2. (A.6) 

An AnB- 



Appendix B: The multifluid approach 

In this section, we compare our results with those obtained by a multifluid 
scheme (arbitrary number of species, but only collisions with neutrals included). 
Neglecting Biermann's battery and diamagnetic terms, Eq. J46) can be written in 
the standard form 



U 
E+ - XB: 



0"H 



1 

lo"ll 



o-p 



(J X B) X B 



B2 



JxB 3_ 

B ^crii 



(B.l) 



where rT|| is the conductivity parallel to the electric field, crp and cth are the 
Pedersen and Hall conductivities, respectively. The conductivities can be written 
in terms of the Hall parameters /J,,, defined as 



qsB \ Ps _ [ ^\ '".; + m„ 



(B.l) 



the ratio of the cyclotron frequency of a particle i and its characteristic frequency 
of collision with neutral particles (notice tha t the Hall parameter is negative for 
negatively-charged particles). From Eq. )46t . the expressions of the conductivi- 
ties result 
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Fig. B.l. Values of t/ad, rjo and?7H as function of the grain radius 
rg for conditions typical of a molecular cloud core (see text). 
The neutral, positive and negative species are H2, HCO^, and 
electrons (long-dashed curves) or grains (short-dashed curves). 
The dash-dotted curves show the same quantities computed with 
the multi-fluid scheme. 
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These expressions are a special case, for two charged fluids, of the general results 
of the multifluid approach. 
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The corresponding resistivities are 
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If the charged particles are electrons and ions, m^ <K m,-. Since the collisional rate 
coefficient ai'e of the same order for electron-neutral and ion-neutral collisions, 
in general \fig„\ 2> /3,„, and one has 



m = LS^Ir/o, 



r]AD = l/3enl/3;n'?0- 



(B.13) 



As an illustration, we analyze the variation of the resistivities appearing in 
Eq. j47t as function of the grain size in a molecular cloud core. 

Consider a molecular cloud core made of neutrals (H2), ions (HCO*), elec- 
trons {e), and negatively -charged dust grains (g"), with radius rg. In particular, 
we assume «„ = 10* cm"', n, = 10"^ cm"', Pg/p = 0.01, a temperature of 
T = 10 K, and a magnetic field B = 0. 1 mG. Each grain is assumed to be spheri- 
cal, with a radius a varying between 10 A and 0. 1 mm, and mean interior density 
2 g cm"' . All grains are assumed to be negatively chai'ged, and the number den- 
sity of free electrons is obtained by the neutrality condition n^ + ng-(rg) = n,. 
For high values of the grain's radius, the number density of grains is small, and 
all the negative charge is canied by electrons. In the opposite limit, the number 
density of small grains is high, a fraction of the grain population is sufficient to 
carry all the negative charge, and the density of free electrons is zero. 
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